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The properties of infinite matter and finite nuclei are studied by using the chiral 
sigma model in the framework of the relativistic mean field theory. We reconstruct an 
extended chiral sigma model in which the omega meson mass is generated dynamically 
by the sigma condensation in the vacuum in the same way as the nucleon mass. All 
the parameters of chiral sigma model are essentially fixed from the hadron properties 
in the free space. In nuclear matter, the saturation property comes out right, but 
the incompressibility is too large and the scalar and vector potentials are about a half 
of the phenomenological ones, respectively. This fact is reflected to the properties of 
finite nuclei. We calculate N = Z even-even mass nuclei between N = 16 and N = 
34. The extended chiral sigma model without the pion mean field leads to the result 
that the magic number appears at N = 18 instead of N = 20 and the magic number 
does not appear at N = 28 due to the above mentioned nuclear matter properties. The 
latter problem, however, could be removed by the introduction of the finite pion mean 
field with the appearance of the magic number at N = 28. We find that the energy 
differences between the spin-orbit partners are reproduced by the finite pion mean field 
which is completely a different mechanism from the standard spin-orbit interaction. 
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1. Introduction 



Chiral symmetry is known to be the most important symmetry in the hadron physics. 
This is because the quantum chromo-dynamics (QCD) is the underlying theory of the strong 
interaction, in which the up and the down quarks have essentially zero masses. Chiral 
symmetry governs the quark dynamics. In the real world, the quarks are confined and chiral 
symmetry is spontaneously broken. As the Nambu-Goldstone boson of the spontaneous 
breaking of chiral symmetry, the pion emerges with almost zero mass. 

At the hadron level, chiral symmetry was described nicely in the linear sigma model 
introduced by Gell-Mann and Levy 1 ). Its non-linear version was proposed by Weinberg 2 \ 
Chiral symmetry and the generation of the hadron mass were described clearly in the Nambu- 
Jona-Lasinio Lagrangian with fermion fields 3 - ) . These Lagrangians have been used for various 
phenomena in hadron physics. We find good description of the pion-nucleon properties in 
terms of these Lagrangians ^ . The pion, which was introduced by Yukawa as the mediator 
of the nucleon-nucleon interaction, received its foundation through the spontaneous chiral 
symmetry breaking 5 \ 

It is then very natural to use the chiral sigma model Lagrangian for the description of 
nuclei. This was performed by several groups in the relativistic mean field approximation 
6 ) " 10 ) . It was found that the use of chiral sigma model in its original form was not satisfactory 
for the description of nuclear matter. An interesting way out was proposed by Boguta et al., 
who introduced a dynamical generation of the omega meson mass in the same way as the 
nucleon mass 7 \ They were able to reproduce the saturation properties of infinite matter with 
the extended chiral sigma (ECS) model. However, the effective mass came out to be large 
and the incompressibility to be very large. The extended chiral sigma model was applied to 
finite nuclei by Savushkin et al. 9 )' 10 ). The binding energies came out to be reasonable, but 
the spin-orbit splitting was too small in the RMF framework. 

Recently, an interesting proposal was made on the role of pion in finite nuclei by some of 
the present authors u \ They made the relativistic mean field calculations with finite pion 
mean field using the TM1 parameter set 12 ). In order to treat the pion mean field, they 
developed the formalism, in which parity mixed single particle states were introduced. With 
the use of the free space pion-nucleon coupling constant, they found that the pion mean 
field becomes finite, especially the effect appears favorably for the jj-closed shell nuclei, and 
the mass dependence of the energy gain associated with the pion behaves as the nuclear 
surface, (V w ) ~ A 2 / 3 . Hence, the name, surface pion condensation, was introduced for this 
phenomenon. 

In this paper, we would like to study the properties of infinite matter in terms of the 
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ECS model by analyzing the non-linear equation of motion for the sigma field and obtain 
the saturation property of nuclear matter. We apply the ECS model to finite nuclei and 
study the properties of the binding energies and the single particle properties. Since the role 
of the finite pion mean field on the binding energies and the spin-orbit splitting has been 
demonstrated in the recent publication u \ we take the formalism to treat the finite pion 
mean field in the RMF framework in the ECS model for the calculation of finite nuclei. We 
would like to study the appearance of the spin-orbit splitting due to the pion mean field by 
studying carefully the single particle spectra of finite nuclei. 

In section 2, we discuss the RMF formalism with the pion mean field. In section 3, 
we study the saturation property of infinite nuclear matter with the original chiral sigma 
model and with the extended chiral sigma model. In section 4, we study finite nuclei with 
the extended chiral sigma model without introducing yet the pion mean field and further 
study the properties of single particle states. We introduce then in section 5 the finite pion 
mean field and discuss the mechanism of the appearance of the magic number effect and 
the energy splittings between the spin-orbit partners. We summarize the present study in 
section 6 together with the statements for the further study. 

§2. Chiral sigma model in the relativistic mean field theory 

We start with the linear sigma model with the omega meson field, which is defined by 
the following Lagrangian 1 ), 

C au = V^d" - g a (a + i l5 r ■ t?) - g^^)V (2-1) 
+ \d»od»o + \d,nd^ - ^(a 2 + Tf 2 ) - + Tf 2 ) 2 

- l -^^ v + ^ V 2 + tT 2 v^ 

+ EG. 

The fields a and tt are the nucleon, sigma and the pion fields, fi and A are the sigma 
model coupling constants. Here we have introduced the explicit chiral symmetry breaking 
term, ea, and in addition the mass generation term for the omega meson due to the sigma 
meson condensation as the case of the nucleon mass in the free space 7 ' . The a — uj coupling 
term of this structure may be derived from the bosonization 13 ) of the Nambu- Jona-Lasinio 
model 3) . 

In a finite nuclear system, it is believed to be essential to use the non-linear representation 
of the chiral symmetry. This is because the pseudoscalar pion-nucleon coupling in the linear 
sigma model makes the coupling of positive and the negative energy states extremely strong 
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and we have to treat the negative energy states very carefully. We can derive the non-linear 
sigma model by introducing new variables and making a suitable transformation, 

a + it ■ 7t = pU, U = e if * ,f " (2-2) 
a + 2 75 r • 7? = pU 5 , U 5 = e^-*//- , 

We further implement the Weinberg transformation for the nucleon field as ip = \fU5\P • We 
obtain then the sigma-omega model Lagrangian in non-linear representation, 

£L = " - 7X - - g ul ^)i> (2-3) 

+ \d» P d»p + ^tr^U^Ut - ^p 2 - ^p 4 

- + l -rJp 2 u^ + epi(U + Ut). 

In the above Lagrangian the vector field, v^ 1 , and the axial vector field, a M , contain the pion 
terms. The vector and the axial vector fields are expanded in terms of the pion field as, 

u* 4 = -^(r ■nr-d^-r -d^nr • vf) + • • •, (2-4) 

8/tt 

The kinetic term is expanded as follows, 

2 2 

^-tr^U^Ut = ^d,nd^ + 0(tt 4 ) + 0(tt 6 ) + • • •, (2-5) 

and the explicitly chiral symmetry breaking term is expanded as follows, 

ep\(V + U f ) = ep(l - 2^vr 2 + ^vf 4 + •■•)• (2-6) 

We take now the lowest order terms in the pion filed and truncate higher order terms. The 
resulting Lagrangian is written as, 

C-L = <M*7^ - 9°P ~ ^-757^ • ~ 9.1^ (2-7) 

- + -gjp 2 u^ + ep(l - ^J2* 2 )- 

We now take the vacuum expectation value for the p field as f w , which is determined by 
the pion decay rate 4 ) , 

(p)o = U (2-8) 
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A new fluctuation field y? may be defined by the equation, 



p = U + tp. (2-9) 
We shall now rewrite the Lagrangian (2.7) in terms of the new field ip, 

= - g a U - 9*<P ~ ttt^t ■ d»n - g^u^ (240) 

+ £(^ + ^)(l--^7T 2 ). 

J TV 

Here, we have dropped a non-essential c-number constant in the above expression. We find 
the term V//tt" is small and drop it as follows, 



1(1 + ^) 2 d^n » l^vf^Tf, (2-11) 

/-, V 9 x 1 E ->2 1 £ ->2 
(1 + — ) 7? PS 7T. 

We have to make the dangerous term, the term linear in tp, zero, which leads to the energy 
minimum condition. 



( e _ mSU)ip — 0, (2-12) 



2 ^ 



Finally the Lagrangian for the new field y? within the above approximations is given as 
follows, 

C' aw = $(i>y^ -M-g a( p- ^757^ ' ~ 9^)i> (2-13) 

1 1 ~ 1 

- -U/u/j^ + -mju^ + gjf^ipuj^ + -gj(p 2 uj ll uj tl , 

where we set M = g a f n , m^ 2 = fi 2 + Xf^ 2 , m 2 = fi 2 + 3A/„- 2 and m u = gZf w . The effective 

mass of the nucleon and omega meson are given by M* = M + g a ip and mj* = + g^(p, 
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respectively. We take the following masses and the pion decay constant as, M = 939 MeV, 
m u = 783 MeV, m n = 139 MeV, and f n = 93 MeV. Then, the other parameters can be 
fixed automatically by the following relations, g a = M/f n = 10.1 and g^, = m^/fn = 8.42. 
The strength of the cubic and quadratic sigma meson self-interactions depends on the sigma 
meson mass through the following relation, A = (m a 2 — m 7r 2 )/2/ 7r 2 , in the chiral sigma 
model. The mass of the sigma meson, m CT , and the coupling constant of omega and nucleon, 
g^, are the free parameters. If we use the KSFR relation for the omega meson 14 ^ ,17 \ and 
the additional relation from the Nambu-Jona-Lasinio model, the mass of the omega meson 
is related to the pion decay constant by m w = (2y / 2/3) f^g^. The factor (2v^2/3) stems from 
the g w = (3/2)</, where g is the universal coupling constant for the vector meson 15 ^' 16 \ As 
we see below, this KSFR relation is very well satisfied in the present model within 6 %. 



§3. Extended chiral sigma model for infinite matter 

We apply first the extended chiral sigma model to infinite matter. It is important to 
reproduce the saturation properties of infinite nuclear matter first. Otherwise, we do not get 
convergence due to the multiple solutions in the Hartree calculation for finite systems. We 
assume that the pion mean field vanishes in infinite matter. Hereafter we write the scalar 
meson field ip in the Lagrangian (2.13) as a, since a is used usually as the scalar meson field 
in the relativistic mean field theory. The equations of motion for the nucleon field and the 
meson fields are written as, 



(i 7/i d" -M-g a a- g u -r°u)ip = 0, 

2 i o \ r 2 i \ 3 ~2 p 2 ~2 2 

m a <j + 3\f w a + Act - g w J^u - g w au = -g a p s , 
mjuj + Ig'Jf^GU + g^cPu = g u p v . 



(3-1) 



with 
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+ M*"-M* z log 



k F + V V + M* 



M* 



(3-2) 



p v = 2k F 3 / (37r 2 ) and the effective mass of the nucleon M* = M+g a a. We note here that now 
the equations of motion of the sigma and omega mesons are coupled due to the dynamical 
mass generation term of the omega meson. This sigma-omega coupling plays an important 
role to obtain reasonable equation of state of nuclear matter. 
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We discuss first the original chiral sigma model for the nuclear matter calculation 7 - 1 . In 
this case, there is no coupling between the equations for the sigma and the omega fields. 
The equation for a is a third order algebraic equation of the sigma together with minus of 
the scalar coupling times the scalar-density, —g a p s , which is a function of the sigma field for 
a fixed density, 

ml<j + 3Xf n a 2 + Xa 3 = -g a p s . (3-3) 

The right hand side increases with decreasing the sigma field and changes sign near a = 
—f n ~ —0.5 fm _1 . We shall focus on the solution above the crossing point, until where the 
effective mass of the nucleon is positive. Below a certain density there appears only one 
solution, while above this density there appear three solutions. We get multiple solutions as 
discussed above. For each solution, there is a corresponding energy, which is not a smooth 
function of the density. Hence, we are not able to get a good behavior for the equation of 
state with the original chiral sigma model. 

The way out to get a good nuclear matter property was suggested by Boguta, who intro- 
duced the dynamical omega meson term 7 -*. The omega mass appears due to the dynamical 
chiral symmetry breaking and hence there is a coupling between the sigma and the omega 
fields. We use this extended chiral sigma model for nuclear matter. The additional term 
provides a pole at the effective nucleon mass being zero, o = —f„, as shown in Fig. 1. Due 
to this reason we find a solution at a small sigma value for each density continuously from 
zero. We are therefore able to obtain a reasonable energy per particle in the entire density 
region for infinite matter. 

In Fig. 2 we provide the energy per particle of nuclear matter as a function of the density 
for the extended chiral sigma model. We take the parameters of the chiral sigma model from 
the properties of mesons as pion mass, m^, omega meson mass, m U) pion decay constant, f n . 
The free parameters, m a and g w , are adjusted to provide the saturation property in the case 
of the extended chiral sigma model. We have fixed the free parameters as, m a = 777 MeV, 
and g u = 7.03. Then, the strength of the cubic and quadratic sigma meson self-interaction 
are fixed as A = 33.8. The saturation properties are the density, p = 0.141 fm -3 , and the 
energy per particle, E/A = -16.1 MeV. We find in this case the incompressibility, K = 650 
MeV. The sigma meson mass chosen here is larger than that used in one boson exchange 
potential, which is around 500 MeV. If we take 500 MeV as the sigma meson mass, the 
attractive force becomes strong and the saturation curve becomes deep. We adjust then 
the omega-nucleon coupling constant, g^, to reproduce the binding energy per particle. The 
energy minimum point appears at quite a small density, p = 0.053 fm~ 3 . The saturation 
condition is not satisfied simultaneously both for the density and binding energy per particle 
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Fig. 1. The equation for a with the a — to coupling term for the case of p = 0.141 fm~ 3 in the 
extended chiral sigma model. There is one solution for each density continuously from the zero 
density. 



using this meson mass. It is interesting to note that the value m a = 777 MeV is very close 
to the one when the chiral mixing angle is chosen at 45° in the generalized chiral model; 
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Fig. 2. The energy per particle of infinite nuclear matter as a function of the density for the 
extended chiral sigma model (solid curve). As a reference the energy per particle in the RMF 
theory with the TM1 parameter set, RMF(TMl), is provided by dashed curve. 
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Fig. 3. The scalar and vector potentials are plotted as a function of the density for the extended 
chiral sigma model shown by solid curve and those for RMF(TMl) by dashed curve. 



As a comparison, the energy per particle of the mean field result with the TM1 parameter 
set is shown together with the present result 12 \ The RMF(TMl) calculation reproduces the 
results of the relativistic Brueckner-Hartree-Fock calculation 19 ) . We see that the present 
equation of state is much harder than the one of RMF(TMl). The incompressibility comes 
out to be 650 MeV, while it is 281 MeV for TM1. In Fig. 3 we plot the vector and the 
scalar potentials and compare with the ones of RMF(TMl). The values are about a half of 
the case of the TM1 parameter set. This is because the extended chiral sigma model has 
solutions at smaller sigma values than those for RMF(TMl). 

We would like to note the consequence of the smaller absolute values of the scalar and 
the vector potentials in finite nuclei as shown in Fig. 3. The summation of the absolute 
values of the scalar and the vector potentials is directly related with the spin-orbit potential 
of finite nuclei. Hence, the fact that these absolute values are about a half of the values of 
RMF(TMl) indicates that the spin-orbit splitting for finite nuclei will come out to be about 
a half of the necessary spin-orbit splittings. 



§4. Extended chiral sigma model for finite nuclei 

We are now in the position to apply the extended chiral sigma (ECS) model, which is 
able to provide the saturation property with the above mentioned features, to finite nuclei. 
For this purpose we take the N = Z even-even mass nuclei to avoid the complication coming 
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from the isovector part of the nucleon-nucleon interaction. We calculate these nuclei using 
the RMF framework with the ECS Lagrangian and compare the results with those of the 
standard RMF calculation with the TM1 parameter set. Since the role of the pion mean 
field on the binding energy and the spin-orbit interaction has been demonstrated in Ref. n ) 
for finite nuclei, we shall introduce the RMF formalism on the treatment of the finite pion 
mean field and study the effect of the finite pion mean field on the nuclear properties. 

We write here the RMF equations for the finite nuclei with the pion mean field. The 
Euler-Lagrange equation gives us the Dirac equation for the nucleon: 

{-ia ■ V + 7o(M + g a a) + gu u + |f 7o7s7 • r°V7r)^ = e^, (4-1) 

and the Klein-Gordon equations for the mesons: 

(-V 2 + m 7r 2 )vr=^p pv , (4-2) 

(-V 2 + m 2 )o = -g aPs - 3Xf n a 2 - Xa 3 + g~J Uoo 2 + <j~W, (4-3) 
(-V 2 + mj)uj = g UJ p v - 2gZ 2 f w (TLd - g'Jo 2 ^, (4-4) 

where we consider the isospin symmetry nucleus, N = Z. There is a symmetry theorem for 
the Hartree-Fock (mean field) approximation with respect to the symmetry of the original 
Lagrangian 20 ^' 21 ). In the isospin symmetric nuclear case, we can verify that the mean field 
Lagrangian is symmetric under the isospin rotation to mix the proton and the neutron states. 
Hence, we can take a special case, where only tt° is finite due to the isospin symmetry of 
the mean field Lagrangian and write it as ir. In fact, we have checked this symmetry by 
performing the mean field calculations with it in one case and with n^ 1 in another case and 
obtained the same energy in both cases 11 ). We take the static approximation and assume 
the time reversal symmetry of the system. We have introduced here qa in the pion nucleon 
coupling in order to fulfill the Goldberger-Treiman relation. In the linear sigma model, we 
get gA — 1. In the mean field approximation, the source terms of the Klein-Gordon equations 
are replaced by their expectation values in the ground state. 

ff V • ^ 7 57tV — > |t-<V • ^757rV) = |f P PV , (4-5) 

^Jtt ^J-k ^Jir 

g a ^ — > 9a(i>^) = gaps, (4-6) 

9w$loil> — > 9u(^hoip) = 9wPv- (4-7) 

The total energy is given by 

£totai = / d 3 rH (4-8) 
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r / 1 , 1 1 a a 

2^ tnjm ~ J d r< -g a p s a + -g^u - -^-rp pv n 

njm \ " J n 



i ^ (o\ -C \ 3 , ^ 4 1~2<- 2 1 ~2 2 2I 

+ g(3A/ w )(T + -o- - -g u f w auj - -g u a u > 
-ZM p -ZM n -E c . m ., 

where we take the center of mass correction as E c . m . = |(4L4.s) MeV. We write here the 
wave functions and the densities for the case of the finite pion mean field. In this case, the 
parity of the nucleon is broken, because the pion source term has the negative parity. The 
nucleon wave functions are then written as, 

^ / iG nKrriT 3^KmC y^T) ^G 'nKm T 'y KmC("*T j 1 ^ 

where the summation over k means the parity mixing, where k is k = —(l^+l) for = j — 1/2 
and k — 1 1 for — j + 1/2. Using these wave functions, we can calculate all the necessary 
densities as, 

Ps ^ ' ^Vnj ~j ^ ] (l^n/tm T I | -^rarem T | \G n Km T \ |-^nfcm T | )? (4-10) 
nj m T 

Pv ^ ' ^Vnj "j ^ ' ( \ GnKrriT \ \Fnnm T \ \Gnnm T \ \FnRm T \ ); (^Tl) 
nj m T 

Ppv = - 2 J2W nj ^-± (4-12) 

nj 



X 



/ A I) 2 } , (G nKmT G nKrriT ) + {G nKmT G nKrriT ) 
[dr r 



~t~ ; (F nKmT F nKrriT ) + (F nKmT F nKUlT 
dr r 



We are now able to calculate the coupled differential equations by doing iterative calculations. 

In this chapter, we discuss first the properties of finite nuclei in terms of the extended 
chiral sigma model without introducing yet the pion mean field. We show the results of 
binding energies per particle of N = Z even-even mass nuclei from N = 16 up to N = 
34 in Fig. 4. We take all the parameters of the extended chiral sigma model as those 
of the nuclear matter (Figs. 2 and 3) except for g u = 7.176 instead of 7.033 for overall 
agreement with the RMF(TMl) results. For comparison, we calculate these nuclei within 
the RMF approximation without pairing nor deformation. The RMF(TMl) provides the 
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20 30 40 50 60 70 80 

A (Mass number) 

Fig. 4. The binding energy per particle for N = Z even-even mass nuclei in the neutron number 
range of N = 16 ~ 34. The binding energies per particle for the case of the extended chiral 
sigma model without and with the pion mean field are shown by the dashed and the solid lines. 
As a comparison, those for the RMF(TMl) are shown by the dotted line. 



magic numbers, which are seen as the binding energy per particle increases at N = Z = 20 
and 28. On the other hand, the extended chiral sigma model without the pion mean field 
provides the magic number behavior only at N = Z = 18 instead of N = Z = 20. 




32 S 36 Ar 40 Ca 44 Ti 48 Cr 52 Fe 56 Ni 60 Zn 64 Ge 68 Se 



Fig. 5. The proton single particle energies for the N = Z even-even mass nuclei in the case of the 
RMF(TMl) theory, where the magic numbers at N = 20 and 28 are visible. 



In order to see why the difference between the two models for the Lagrangian arises, 
we show in Fig. 5 the single particle levels for the two models. In the case of the TM1 
parameter set shown in Fig. 5, the shell gaps are clearly visible at N = 20 and 28. The 
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Fig. 6. The proton single particle energies for the N = Z even-even mass nuclei in the case of the 
extended chiral sigma model without the pion mean field. lsi/2 orbit is pushed up and the N 
= 20 magic number is shifted to N = 18. The spin orbit splitting between 0/ 7 / 2 and 0/ 5 / 2 is 
small and the magic number at N = 28 is not visible. 



magic number at N = 20 is due to the central potential, while the magic numbers at N = 28 
comes from the spin-orbit splitting of the Of-orbit. This is definitely due to the fact that the 
vector potential and the scalar potential in nuclear matter are large so as to provide the large 
spin-orbit splitting. On the other hand, the single particle spectrum of the extended chiral 
sigma model is quite different from this case as seen in Fig. 6. Most remarkable structure is 
that the lsi/2 orbit is strongly pushed up. Due to this reason the 0^3/2 orbit becomes the 
magic shell at N = 18 and the magic number appears at N = 18 instead of N = 20. We see 
also not strong spin-orbit splitting and hence there appears no shell gap at N = 28. The 
first discrepancy could be due to the large incompressibility as seen in the nuclear matter 
energy per particle as seen in Fig. 2. The other is due to the relatively small vector and 
scalar potentials in nuclear matter as seen in Fig. 3. 

We would like to detail further the discussion on the spin-orbit splitting in the situation 
where the compressibility is very large as the ECS model, since the spin-orbit splitting is 
related with the behavior of the scalar-vector potential difference in the surface region. We 
show in Fig. 7 how the scalar-vector potential difference behaves as a function of r, which 
is defined as U is — U v — U S - The magnitude of the ECS model is about a half of the TM1 
case. The spin-orbit potential is then denned by eliminating the small component in the 
relativistic wave function and by getting the spin-orbit operator explicitly as, 

K -f(«H-f-W r? <4 ' 13) 
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Fig. 7. The scalar- vector potential difference, Ui s , which is related with the spin-orbit potential, 
as a function of the radial coordinate, r. The potential difference for the case of the TM1 
parameter set is shown by the dashed curve and for the ECS model is shown by the solid curve. 

The spin-orbit potential is proportional to the derivative of the scalar-vector potential dif- 
ference, which emphasizes the contribution from the nuclear surface. Hence, to compare the 
magnitude of the spin-orbit effects of the two cases, we calculate the volume integrals of Vi s ; 



We use for the e the value corresponding to the binding energy of 8MeV, and obtain the 
ratio of the two cases as 0.48, which is again about a half. Hence, the spin-orbit effect for 
the ECS model is about a half of the TM1 case, which could be seen already in the single 
particle spectra shown in Fig. 6. 



We include now the pion mean field in the relativistic mean field calculation 22 )- 23 ). The 
method of the numerical calculation is provided in the paper of Toki et al. n ) and Sugimoto 
et al. 18 ). The results on the binding energy per particle are shown in Fig. 4. In this 
calculation we take 1.15 instead of the experimental axial coupling constant g\ = 1.25 due 
to the Goldberger-Treimann relation. We take this smaller value in order to reproduce the 
binding energy for 56 Ni. It is very interesting to see that the magic number effect at N = 28 
appears as the binding energy per particle increases at N = 28. This large effect of the finite 
pion mean field for the jj-closed shell nuclei has been demonstrated in the previous work 11 ). 




(4-14) 



5. Finite pion mean field for finite nuclei 
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Fig. 8. The proton single particle energies for the N = Z even-even mass nuclei in the case of 
the extended chiral sigma model with the pion mean field. The spin-orbit splitting is made 
large due to the finite pion mean field, which is visible as centered at the N = Z = 28 nucleus. 
We note that while the total angular momentum is a good quantum number, but the angular 
momentum is not exact, we write the dominant angular momentum beside each single particle 
state. 




Fig. 9. The schematic picture of the single particle states and the occupied particles in the 56 Ni 
nucleus. 



We give here an intuitive explanation to understand the energy curve of the magic struc- 
ture in Fig. 4 to be given by the finite pion mean field using the schematic picture in Fig. 
9. To proceed, we have to know first the effect of the finite pion mean field in terms of the 
shell model. The discussion of the parity projection, done in the previous publication n \ 
clearly shows that the pionic correlations due to the finite pion mean field is expressed by 
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the coherent 0~ particle-hole excitations, in which the coupling of the different parity states 
I and V — I ± 1 with the same total spin j in the shell model language. In the discussion 
of the contribution to the pionic correlations from various single particle states, the highest 
spin state in each major shell has a special role. Only this highest spin state does not find 
the partner to form the 0~ state in the lower major shells. However, if this state is filled 
by nucleons, those nucleons are able to find the 0~ partners in the higher major shells by 
making particle-hole excitations. Hence, the position of the highest spin state in a major 
shell with respect to the Fermi surface is important for the strength of pionic correlations in 
nuclei. 

In the case of discussion, the highest spin state is the /t/ 2 state as shown in Fig. 9. In 
the 40 Ca case, the occupied states can not couple with the state to form 0~ and the f 7 / 2 
level is not used at all for the pionic correlations. In the next 44 Ti case, nucleons start to 
occupy the f?/2 level, and these nucleons are used for the 0~ particle-hole excitations into 
97/2 levels. The number of particles to be used by the pionic correlation increases as the 
nucleon number is increased until 56 Ni, where the f-j/2 level is completely occupied. For the 
nuclei above 56 Ni, the upper shells as / 5 / 2 are to be occupied and those states are not used 
for 0~ particle-hole excitations from the d$/2 level below caused by the pionic correlation 
due to the Pauli blocking. For 56 Ni, the pionic correlation becomes maximum. This is the 
reason why 56 Ni obtains the largest pionic correlation energy, which leads to the appearance 
of the magic number at N = 28. 

We discuss now the effect of the finite pion mean field on the single particle energies. 
We show in Fig. 8 the single particle spectra for various nuclei. We see clearly the large 
energy differences between the spin-orbit partners to be produced by the finite pion mean 
field as the energy differences become maximum for nuclei at N = 28. The pion mean field 
makes coupling of different parity states with the same total spin. The Osi/ 2 state repels each 
other with the 0pi/ 2 state and therefore the 0si/ 2 state is pushed down and the 0pi/ 2 state 
is pushed up. The next partner is 0p 3 / 2 and 0d 3 / 2 . The 0p 3 / 2 state is pushed down, while 
the 0d 3 / 2 state is pushed up. The next partner is 0d 5 / 2 an d 0/ 5 / 2 . The 0d 5 / 2 state is pushed 
down, while the 0/ 5 / 2 state is pushed up. This pion mean field effect continues to higher spin 
partners. This coupling of the different parity states with the same total spin due to the 
finite pion mean field causes the splittings of the spin-orbit partners as seen clearly for the 
Op spin-orbit partner, Od spin-orbit partner and 0/ spin-orbit partner in 56 Ni. It is extremely 
interesting to see that the appearance of the energy splitting between the spin-orbit partners 
for the case of the finite pion mean field is caused by completely a different mechanism from 
the case of the spin-orbit interaction. 

We would like to see the contributions of each term in the Lagrangian for the cases with 
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and without the pion mean field in Table 1. 



Table I. The binding energy per particle (BE/ A) and the contributions of the sum of sigma and 
omega, (U CT + U w ), kinetic (KE), pion (U^), non-linear term (NL), sigma-omega coupling term 
(CP) and Coulomb (Uc) energies per nucleon in MeV for 56 Ni in the extended chiral sigma 
model. 





BE/ A 




KE 


u* 


NL 


CP 


U c 


with 7r field 


8.6 


-21.8 


20.9 


-2.9 


8.1 


-15.4 


2.6 


without 7r field 


8.4 


-22.6 


18.8 





8.0 


-15.0 


2.6 



The binding energy increases slightly by making the pion mean field finite. The pion term 
contributes attractively and the energy gain due to the pion term is obtained by making the 
kinetic energy and the sum of the sigma and omega potential terms increase. The structure 
of the wave functions changes largely, while the total energy is kept almost unchanged. This 
change of the structure will make the observables associated with the spin quantities change 
largely. The effect of the structure change on various observables will be studied in the near 
future. 



§6. Conclusion 

We have studied infinite nuclear matter and finite nuclei with the nucleon number N = 
Z even-even mass in the range of N = 16 and N = 34 using the chiral sigma model, which 
is good for hadron physics. The direct application of the chiral sigma model is not able to 
provide the good saturation property of infinite matter. We have then used the extended 
chiral sigma (ECS) model, in which the omega meson mass is dynamically generated by 
the sigma condensation as the nucleon mass. This ECS model is able to provide a good 
saturation property, although the incompressibility comes out to be too large. Another 
characteristic property of the ECS model is that the scalar and vector potentials are about 
a half of the case of the RMF(TMl) model in nuclear matter. 

We have then applied this ECS model to finite nuclei. The ECS model without the pion 
mean field gives the result that the magic number appears at N = 18 not at N = 20. This 
result comes from the large incompressibility found in the equation of state as K = 650 
MeV. This property of the ECS model provides the mean field central potential repulsive in 
the interior region and the ls-orbit is extremely pushed up. Due to this, the magic number 
appears at N = 18 instead of N = 20. We note that this problem originates from the ECS 
model treated in the present framework and the finite pion mean field under the mean field 
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approximation does not remove this difficulty. There are several possibilities to be worked 
out to cure this problem as the effect of Dirac sea, the parity projection, and the Fock term. 

The ECS model without the pion mean field provides the result that the magic number 
does not appear at N = 28. This result comes from another characteristic property of the 
ECS model, which is the small scalar and vector potentials in nuclear matter. The scalar and 
vector potentials lead directly to the strength of the spin-orbit interaction in finite system. 
Since the spin-orbit interaction given by the ECS model is about a half of those of the 
standard RMF calculation with the TM1 parameter set, the energy splittings between the 
spin-orbit partners are small and, therefore, there appears no magic effect at N = 28. As 
for this point, it is important to introduce the pion mean field by breaking the parity of 
the single particle states in the ECS model Lagrangian. Since the role of the pion mean 
field on the jj-closed shell nuclei has been demonstrated in the previous publication n \ we 
have introduced the parity mixed intrinsic single particle states in order to treat the pion 
mean field in finite nuclei. We followed the formulation of Sugimoto et al. 18 - ) in the RMF 
framework. We have found that the magic number effect appears at N = 28. We have 
studied the change of the single particle spectrum due to the finite pion mean field. It is 
extremely interesting to find that the spin-orbit partners are split largely by the pion mean 
field effect. Namely, the parity partners as (si/2 andpi/2), G°3/2 and ^3/2) and (^5/2 and $5/2) 
are pushed out each other due to the pion mean field and as the consequence the spin-orbit 
partners are split largely like the ones of the ordinary spin-orbit splittings. This is related 
with the energy differences of the spin-orbit partners caused by the energy loss of the tensor 
(pionic) correlations due to the Pauli blocking 24 \ 

It is gratifying to observe that first the extended chiral sigma model, which has the chiral 
symmetry and its dynamical symmetry breaking, is able to provide the nuclear property 
with only a small adjustment of the parameters in the Lagrangian. The energy splitting 
between the spin-orbit partners appears remarkably in the ECS model with the pion mean 
field. The most important consequence obtained in this study is that this energy splitting 
is caused by the pion mean field which is completely a different mechanism from the case 
of the spin-orbit interaction introduced phenomenologically. This suggests the origin of the 
magic effect of jj-closed shell nuclei. 

Acknowledgement 

We acknowledge fruitful discussions with Prof. Y. Akaishi, Prof. H. Horiuchi and Prof. 
I. Tanihata on the roll of the pion in nuclear physics. We are grateful to Dr. D. Jido 
for helpful discussions on the linear sigma model. We thank Prof. E. Oset for reading 



18 



manuscript and valuable discussions. This work is supported in part by the Grant-in-aid for 
Scientific Research (B) 14340076 of the Ministry of Education, Culture, Sports, Science and 
Technology of Japan. 



References 

M. Gell-Mann and M. Levy, Nuovo Cimento 16 (1960), 705. 

S. Weinberg, Phys. Rev. 166 (1968), 1568; 177 (1969), 2604. 

Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961), 345; 124 (1961), 246. 

B. W. Lee, Chiral Dynamics, Gordon and Breach Science publishers. 

H. Yukawa, Proc. Phys.-Math. Soc. Jpn. 17 (1935), 48. 

J.D. Walecka, Ann. of Phys. 83 (1974), 491; B.D. Serot and J.D. Walecka, in Ad- 
vances in Nuclear Physics, edited by J.W. Negele and E. Vogt (Plenum Press, New 
York, 1986), vol. 16, p. 1. 

J. Boguta, Phys. Lett. 120B (1983), 34; 128B (1983), 19. 

J. Kunz, D. Masak, U. Post, and J. Boguta, Phys. Lett. 169B (1986), 133. 

V. N. Fomenko, S. Marcos, and L.N. Savushkin,J. of Phys. G19 (1993), 545. 

V.N. Fomenko, L.N. Savushkin, S. Marcos, R. Niembro, and M.L. Quelle, J. of Phys. 

G21 (1995), 53. 

H. Toki, S. Sugimoto and K. Ikeda, Prog. Theor. Phys. 108 (2002), 903. 
Y. Sugahara and H. Toki, Nucl. Phys. A579 (1994), 557. 
A. Hosaka et al., to be published (2003). 

K. Kawarabayashi and M. Suzuki, Phys. Rev. Lett. 16 (1996), 255. 
A. Hosaka, Phys. Lett. 244B (1990), 363. 

O. Kaymakcalan, S. Rajeev, and J. Schechter, Phys. Rev. D30 (1984), 594. 

Riazuddin and Fayyazuddin, Phys. Rev. 147 (1966), 1071. 

S. Sugimoto, H. Toki, and K. Ikeda, to be published. 

R. Brockmann and R. Machleidt, Phys. Rev. C42 (1990), 1965. 

G. Ripka, Adv. Nucl. Phys. 1 (1968), 183. 

H. Horiuchi and K. Ikeda, Int. Rev. Nucl. Phys. 4 (1985), 1. 
H. Toki and W. Weise, Phys. Rev. Lett. 42 (1979), 1034. 

E. Oset, H. Toki, and W. Weise, Phys. Rep. 83 (1982), 281. 

S. Takagi, W. Watari, and M. Yasuno, Prog. Theor. Phys. 22 (1959), 549; T. Tera- 
sawa, Prog. Theor. Phys. 23 (1960), 87; A. Arima and T. Terasawa, Prog. Theor. 
Phys. 23 (1960), 115. 



19 



